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We develop a detailed microscopic theory describing dc Josephson effect and Andreev bound 
states in superconducting junctions with a half-metal. In such systems the supercurrent is caused by 
triplet pairing states emerging due to spin-flip scattering at the interfaces between superconducting 
electrodes and the half-metal. For sufficiently clean metals we provide a detailed non-perturbative 
description of the Josephson current at arbitrary transmissions and spin-flip scattering parameters 
for both interfaces. Our analysis demonstrates that the behavior of both the Josephson current and 
Andreev bound states crucially depends on the strength of spin-flip scattering showing a rich variety 
of features which can be tested in future experiments. 
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I. INTRODUCTION 

Recent experiments^ strongly indicate the possibility 
to realize non-vanishing supercurrent across sufficiently 
thick half-metal (H) layers embedded in-between two s- 
wave BCS superconductors (S). This physical situation 
appears rather non-trivial. Indeed, in conventional SNS 
junctions (N stands for spin-isotropic normal metal) the 
supercurrent is carried by (spin-singict) Cooper pairs 
penetrating into the N-metal layer from both supercon- 
ductors due to the proximity effect. However, half- metals 
are fully spin polarized materials acting as insulators for 
electrons with one of the two spin directions. Hence, 
penetration of spin-singlet electron pairs into half met- 
als should be prohibited and no supercurrent would be 
possible. 

Recently it was realized^ that this situation changes 
qualitatively if one allows for spin-flip scattering at HS 
interfaces. Such scattering enables conversion of spin- 
singlet pairing in S-electrodes into spin-triplet pairing in 
a half-mctal. In this way superconducting correlations 
can survive even in a half- metal ferromagnct, thus "un- 
blocking" the supercurrent across SHS junctions. Sub- 
sequent numerical and analytical studies^i^i^ confirmed 
this physical picture also extending it to structures with 
disorder. In particular, it was argued^ that depending 
on the degree of disorder in the H-metal the origin of 
triplet pairing there can change from the p-wave type 
to the odd-frequency one^. It was also realized that the 
presence of triplet pairing in strong ferromagnets and in 
half-metals can cause the so-called 7r-junction behavior 
of the system!. 

Despite all these important developments the issue is 
yet far from settled. The goal of this paper is to work out 
a complete theory of dc Josephson effect in clean SHS het- 
erostructures at arbitrary transmissions of HS interfaces. 
Electron scattering at these interfaces will be described 
by the most general scattering matrices which fully ac- 
count for spin-flip processes. Depending on the system 



parameters we will find a rich variety of the results both 
for the temperature dependence of the supercurrent and 
for the current-phase relation in SHS junctions demon- 
strating crucial importance of spin-flip scattering at HS 
interfaces. In addition we will also address Andreev level 
quantization and show that this phenomenon in SHS 
junctions acquires qualitatively new features which are 
not present in conventional SNS structures. 

The structure of the paper is as follows. In Sec. [U we 
will employ the quasiclassical formalism which enables us 
to exactly evaluate the Josephson current in SHS struc- 
tures with many conducting channels at any transmis- 
sions of SH interfaces. In Sec. IIIII we will develop a 
more general approach which also accounts for resonant 
effects and allows to include structures with few conduct- 
ing channels into consideration. Within this approach we 
will describe both the Josephson current and Andreev 
levels in SHS junctions and establish the correspondence 
to the results derived in Sec. |TT1 In Sec. llVl we will briefly 
summarize our main observations. 



II. QUASICLASSICAL ANALYSIS 



In this section we will consider a general model of a 
clean SNS junction with spin-active interfaces. We will 
then specify the scattering matrices of NS interfaces in a 
way appropriate to describe SHS heterostructures. Here 
we will treat the systems with many conducting chan- 
nels. For this reason it will be sufficient to employ the 
quasiclassical formalism of energy-integrated Matsubara 
Green functions^i^. 
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FIG. 1: SNS junction and Riccati amplitudes in the clean limit. The functions 7^, Vi 7; Ti are Riccati amplitudes at the 
corresponding NS interface. 7± and 7± are Riccati amplitudes in the middle of the normal metal layer. Quasiparticle momentum 
directions are indicated by arrows. 



A. Riccati parameterization 

In the ballistic limit the Eilenberger equations take the 
form 

iwnfa - A(r), 'g{pF, w„, r) + ivp'^giPF, <^n, r) = 0, 

where [a,b] = db — ba, u>n = ■nT(2n + 1) is Matsubara 
frequency, Pf = rnvF is the electron Fermi momentum 
vector and is the Pauli matrix in Nambu space. The 
function g also obeys the normalization condition g^ = 1. 
Green function g and A are 4x4 matrices in Nambu and 
spin spaces. In Nambu space they can be parameterized 
as 



9 f 
f ~9 



A 



Aia2 
A*ic72 



(2) 



where g, /, /, g are 2x2 matrices in the spin space, A is 
the BCS order parameter and at are Pauli matrices. For 
simplicity we will only consider the case of spin-singlet 
isotropic pairing in superconducting electrodes. As usu- 
ally, the superconducting order parameter in the normal 
layer is set to be equal to zero. 

The equilibrium current density is defined by the stan- 
dard relation 

j(r) = eNonT ^ Im (t;^ Sp[f3g(pi., r)]) , (3) 

UJ„>0 

where A'o = TOpi?/27r^ is the density of states at the Fermi 
level and angular brackets (...) denote averaging over the 
Fermi momentum. 

The above matrix Green functions can be conveniently 
parameterized^" by the two Riccati amplitudes 7 and 7: 



(l-77)-i(l+77) 
-2(l-77)-i7 



2(1"77)"'7 
-(l-77)-i(l+77) 



(4) 



the following set of decoupled equations for Riccati 
amplitudes'*^ 

Wi?V7 + 2iciJ7 ~ 7A*icr27 — Aicr2, (5) 
Wi?V7 — 2iciJ7 = 7Ai(T27 — A*icr2- (6) 

Solving Eqs. ([5]) and ^ inside the normal metal we 
obtain following relations between Riccati amplitudes 

7+ = Fi/ exp{-LL>nd/\vFx\) = 72' ex.p{uJnd/\vFx\), (7) 

7_ = 71, ex.p{ujnd/\vFx\) = r2' ex.p{-uj„d/\vFx\), (8) 

7+ = 71' exp{ujnd/\vFx\) = f 2' ex.p{-uJnd/\vFx\), (9) 

7- = f 1/ ej<.p{-uJnd/\vFx\) = 72' e:Kp{ujnd/\vFx\), (10) 

where 7^' , Ti> , 7^/ , F^/ are Riccati amplitudes at the corre- 
sponding interface and 7± and 7± are Riccati amplitudes 
in the middle of the normal metal slab (see Fig. [T]for de- 
tails), d is a distance between two NS interfaces. 

Deep inside the superconducting electrodes we apply 
the following asymptotic conditions 

71 ^ -c^2a(w)e"*x/^ 71 = o•2a(w)e''^/^ (11) 

72 = -a2a{uj)e'^/^ 72 = <J2a{uj)e-'^/^ , (12) 



With the aid of the above parameterization one can iden- 
tically transform the quasiclassical equations ([1]) into 



where a{uj) — uP- + |Ap — ti;)/|A| and x is the su- 
perconducting phase difference across the junction. Here 
and below we assume that the order parameter is spa- 
tially uniform inside superconducting electrodes. This 
choice can always be parametrically justified by assum- 
ing the proper junction geometry and/or interface trans- 
mission values. For simplicity we also assume that the 
absolute values of the superconducting order parameter 
are identical in both S-electrodes. Generalization of our 
approach to the case of anisotropic pairing and asymmet- 
ric electrodes with |Ai| 7^ IA2I is straightforward. 



B. Boundary conditions at metallic interfaces 

The above quasiclassical equations should be supple- 
mented by appropriate boundary conditions at the inter- 
faces. In the case of specularly reflecting spin-degenerate 
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interfaces these conditions were derived by Zaitse^iii and 
later generalized to spin- active interfaces in Ref. [H- 

Similarly to Ref. [Ij it will be convenient for us to use 
the boundary conditions formulated directly in terms of 
Riccati amplitudes. Let us consider the first NS interface 
and explicitly specify the relations between Riccati am- 
plitudes for incoming and outgoing electron trajectories, 
see Fig. [TJ For instance, the boundary conditions for Fi, 
can be written in the forn>iiii^ 



where we introduced the following 2x2 matrices 



Fi = riaiS_ii + tir/vS_lj^,, 
f 1 = i^nTl^lr + S;^i7i'^ir- 



(13) 
(14) 



Here we defined the transmission (t) and reflection (r) 
amplitudes as: 



where 



ru = [(3^,lS+-p-,l,S+,]-'Pvl 



(15) 
(16) 
(17) 
(18) 

(19) 



Matrices Sij and 5^ are building blocks of the full elec- 
tron and hole interface S-matrices^^ 



S 



Sii Siy 
Svi Syy 



s 



V^i'i ^I'l' 



(20) 



Boundary conditions for Fi/, Fi' can be obtained from 
the above equations simply by replacing 1 1'. Bound- 
ary conditions describing electron scattering at the sec- 
ond interface are formulated analogously. 

Combining the above relations between Riccati ampli- 
tudes we obtain matrix quadratic equations for the ma- 
trices 7]± = —i^±/j2 and f)± = i<72"f± 

?7+|c2ai exp{-ujnd/\vFx\) + a2bi exp{uJnd/\vFx\)^r]+~ 
- r;+|c2Ci exp{-2uJnd/\vFx\) + a2rfi| + 
-I- |d2ai + 62^1 exp{2ujnd/\vFx\)^ri+- 
- |d2Ci exp{-ujrid/\vFx\) + b2di exp(w„d/|i;Fj;|)| = 0, 



(21) 



77+ 1 d2 ci exp {-uj„d/\vFx\) + b2di exp {ujnd/\vFx\)^n 
- ry+|(i2ai + fe2&i exp{2uJnd/\vFx\)^ + 
+ |c2Ci exp{~2und/\vFx\) + a2rfi|?y+- 
- |c2ai exp{-LUnd/\vFx\) + a2bi exp{ujnd/\vFx\)^ 



= 0, 
(22) 



ai = [i(y2^iiniPii ^i'i\ ' 




(23) 


Ol - [ifyy - Oll/Pll OpiJ 




(24j 


ci — [<^2{3.vv +j2ii'7iPii li^v 




(25) 


di = [Jn'Pii 7ii2.i'i«<^2j , 




(26) 


02 = [icr2 5^2' 72/3^2^ 5'^;2] , 




(27) 


62 — \S2'2' ~ 'S'22'/'22 '^2'2] 




(28) 


C2 = [(^2{S+,^, + S+^,l2(i22l2Sl 


'2)^2] , 


(29) 


d2 = [S22'P22^l2S_^,2'^'^2] ■ 




(30) 



Matrix equations for the 77^ and 7)_ can be obtained 
from Eqs. (|2T|) and ((22)) by substituting 77+ — > 77_, 77+ ^ 
7y_ and interchanging of indices 1 <-!■ 2. 

For arbitrary interface S'-matrices the matrix equations 
pT|) and ((22)) can be reduced to scalar quartic equations 
with very cumbersome general solutions. These solu- 
tions, however, become simpler for some particular in- 
terface models. Here we will stick to the case of SHS 
junctions in which wc should specify the scattering S*- 
matrix for the interface between spin-isotropic normal 
metal (or superconductor) and fully spin polarized ferro- 
magnet. We will demonstrate that in this case Eqs. ((^T)) 
and ((22)) can be solved in a transparent and compact way. 



C. Scattering matrices 

Let specify the scattering S'-matrices for both SH in- 
terfaces. For simplicity we will assume that these matri- 
ces depend only on the incidence angle but not on the 
azimuthal one. Then the following relation between elec- 
tron and hole 5-niatrices holds: S — . For a half- 
metal slab between two superconducting electrodes the 
corresponding interface scattering matrices contain 3x3 
nontrivial sub-matrices, i.e. 



S 








(31) 



It is straightforward to demonstrate that the Joseph- 
son current is invariant under the following transforma- 
tion of the S'-matrices: 



Si 

52 

5, - 



Ui 
V 



Si 





v+ 



u*i A /c/f 
v*)-^\ Q v'^ 



U2 







V v+ 
V* -A 



(32) 
(33) 
(34) 
(35) 
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where Ui,U2:V e SU{2). With the aid of the above 
transformation we can always reduce the iS-matrices to 
the following form 



S = 



V 



ae*^ cos ly 



ae 



-/3*e*'' COS a*e*^ 
— e*''sini' 




-f3*e- 



sin v 




e cos ly 




(36) 



V 



where |ap + = 1 and ly, ( are real. We also note 
that in general the first and the second interfaces are 
characterized by two different sets of parameters a, /?, 
v and i.e. these interfaces are described by different 
scattering matrices. 

Here and below the parameter sin^ v defines an effec- 
tive interface transmission. The normal state differential 
conductance (dl/dV) for the metallic interface described 



by the 5-matrix (|36p is proportional to sin 



Gnn — 



2tt 



Ipii I<pf 



■ 2 

sm ly 



27r ^ 

k 



^ sin^ ly (37) 



where the index k labels conducting channels of our junc- 
tion. The limit v — corresponds to completely impen- 
etrable interfaces. The parameter P is responsible for 
spin-flip scattering of electrons at the interface. 

Finally, we point out that by virtue of Eqs. (|32|) - (|35l) 
the scattering matrix employed in the analysis of Ref. 
can be reduced to our expression ([36)1 provided we iden- 
tify 



cos !/ = f — 



2W 



sm v 



t 



TT 



j-2 I 4-2 



je/2 _ 



j-2 I j-2 

^ XT 



/3 = 21^11^ sin(0/2)e-*(«iT+fn), ^ ^ q, 
tff -l- 1 1 1- 



(38) 
(39) 
(40) 



where we use the notations from Ref. 



D. Josephson current 



Now we are ready to evaluate the Josephson current. 
In SNS systems with interface 5'-matrices of the form 
(I36p the matrices ai^hi,Ci,di,'q±,f}± have the following 
structure 





01 



d, 




h 

1 

77± 




»7± 



I 

c, 


fl± 



(41) 
(42) 



Since all the matrices ai,bi,Ci,di,r]±,fi± have only one 
nontrivial matrix element it suffices to denote this ele- 
ment by the same symbol as the corresponding matrix 
itself. Then from Eqs. (PT|) . ([221) we derive a simple 
quadratic equations for the scalar variables ry± and r/i- 
Resolving this equation and making use of the Riccati 
parameterization described in Sec. Ill Al and lllBl wc con- 
struct the Green-Eilenberger function for our junction. 
Substituting this function into Eq. Q we arrive at the 
following general result for the Josephson current 



I{X) = -^eAT J2 



aj„>0 



(2^Y 



P\\\<PF 



.2.2 . . ^^^^ 

where I?i2 — |/3i||/32| sin j^i sin V2, A is the junction 
cross section, 



Q{uj) = 



2a^(l — a^)^|/3i||/32| sin^ vi siv? V2 cosx— 
- ((1 - a^){\ - cos^ vx) -f c? \axf) ((1 - a^)(l - cos^ V2) + l^al^) exp(2w„d/|t;i.x|) 
((1 - a^){-a^ + cos^ vr) ^ \axf^ ((1 - c?){-a^ + cos^ V2) -f |a2|^) exp(-2w„d/|wF:E|) 



1/2 



4 1(1 - a^)^ cosi^i a^ail \(\ - c?f C0&V2 ^ c? \ , (44) 



and We observe that the Josephson current is proportional to 

the parameter I?i2 which contains effective transmissions 

X = X + 3C2 + arg/32 -3Ci -arg^i, = ai + a* cosu^. 

(45) 
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of both interfaces sin^ vi sin^ 1/2 and the two spin-flip fac- 
tors \ and 1/32 1- The scattering matrix parameters Cij 
arg /?! , C2 , and arg (32 enter into our result only in com- 
bination with the phase difference x PSj) . 



Let us briefly analyze the above general expression for 
the Josephson current. At small transmissions i^i, 1^2 ^ 
and for sufficiently long junctions d S.oi'^i + ^2) 
Eq. (gS]) reduces to 



eAT ^ 

UJ„>0: 



d'PW 



bii l<P-F 



i2n) 



A^iul\(3,\\P2WMsinx 



A2(Reai)2][cj2 + A2(Rea2)2] smh{2ujnd/\vF^ 



(46) 



which matches with the analogous result derived in Ref. 
provided we set Reai,2 — 1- For vp/d <C T <C A the 
integral in Eq. (|43|) is dominated by the contribution 
of momenta values sufficiently close to p\\ = and the 
dependence for the Josephson current on the junction 
thickness d acquires the standard exponential form 



1 — TTT-, exp(— zttJ d/VF) — — 2~Z — 2 



ai \a2\ 



irA^d 

P||=0 

(47) 

At lower temperatures T <C w^p/d < A Eq. ([33]) yields 
the power-law dependence on d: 



I : 



7C{3)eA 



47rd3A2 J (27r)2 

|P|i l<PF 



lail^ |a2p 



(48) 



I.e. 



/ cx l/d"^ at T 0. In the limit of the small trans- 
parencies ^1^2 ^ 1 or spin-flip factors |/32| ^ 1 the 
term in Q{uj) proportional to cosx is irrelevant, and the 
current-phase relation becomes purely sinusoidal 



Hx) = -IcSin{x- Xo), 



(49) 



where the phase shift xo depends on the scattering ma- 
trix parameters (j45p and, hence, in general can take any 
value. Provided these values change randomly for differ- 
ent conducting channels the net Josephson current across 
the system can be significantly reduced. For symmetric 
interfaces the phase shift xo is identically zero and the 
TT-junction behavior is realized. 

Note that the expression (|46p formally diverges in the 
limit of small d illustrating insufficiency of the perturba- 
tive (in the transmission) approach for the case of suf- 
ficiently short SHS junctions. This divergence is, how- 
ever, regularized within the non-perturbative approach 
adopted here. From Eqs. (|43)) . ((44l) we observe that for 
very short junctions d ^ Co('^i + 1^2) (and for Reai — 
Re 02) the Josephson current scales with the tunnel in- 
terface transmissions as 



/cx 



(50) 



Obviously this dependence cannot be derived within a 
simple perturbative approach in i'i,2- 




-0.012 



-0.014 



sin fi = sin 1^2 = 0.5 
T = OmTa, d = 5{o 



Al=/32| = 1 - 
ftl = IftI = 0.99- 
ftl = lAI = 0.95- 



0.5 



1.5 



2.5 



FIG. 2: Phase dependence of the Josephson current at low 
temperatures and different spin-flip factors /3. For simplicity 
we consider the case of identical interfaces and assume that 
scattering parameters v, /3, a, C, are momentum independent. 
Parameter a is chosen to be real. The phase coherence length 
^0 is defined as ^0 = vf / {2nTc). 



A detailed analysis of analytical expressions for the 
current in the limit of high interface transmissions will be 
postponed to the next section. Here we only present Figs. 
[2] and [3] illustrating some key features of our general re- 
sults (H51) - (|45l) . The current-phase relation for sufficiently 
thick SHS junction is depicted in Fig. H at T = O.OITc 
for the case of highly transmitting identical interfaces. 
We observe that within the interval of Josephson phases 
X ranging from zero to tt the current is always negative. 
Strong deviations from the sinusoidal current-phase rela- 
tion emerge only provided the spin-flip parameter |/3| is 
very close to unity which corresponds to (almost) com- 
plete spin-flip scattering at both interfaces. Temperature 
dependence of the critical Josephson current is shown in 
Fig. [3] for different values of the spin-flip parameter |/3|. 
Similarly to Refs. ii in a wide parameter range we find 
non-monotonous dependence of the critical current on T 
with a maximum typically below 0.2Tc. We also observe 
that this feature disappears as the spin-flip factor \(3\ ap- 
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0.014 



l/3i|=A| = l 

|/3i| = lAI = 0.99 

|/3i| = lAI = 0.95 

sin^ vi = sin^ U2 = 0.5 

d = 5{o 




0.002 



0.4 0.6 

T/Te 

FIG. 3: Critical Josephson current Ic as a function of tem- 
perature for different spin-flip factors /3. The inset shows the 
current-phase relation for fHi = (32 = 0.99 and at temperatures 
T/Tc^O.Ol (solid line), 0.1 (dashed line), 0.2 (dotted line). 
For simplicity we consider the case of identical interfaces and 
assume that scattering parameters v, (3, a, are momentum 
independent. Parameter a is chosen to be real. 



A. General formalism 

Below we will make use of the more general micro- 
scopic formalism based of the Gorkov equations. Let us 
introduce the matrix Matsubara Green functions"'^* 

Gw{xi,X2,Ti - T2) = - {Tr1pl{xi,Ti)iji,{x2,T2)) , (51) 

where the indices are numbered as 

tpl=1ph i^2=1pl, Ips^lpp TpA^tpp (52) 

As before, assuming singlet pairing in superconducting 
electrodes, one can write down the standard Gorkov 
equations 

[iun - [e(a;i) + A(a;i)]) G{xi,X2,ujn) = 5{xi - X2), 

. 

where we performed the Fourier transformation with re- 
spect to Ti — T2 introducing the Matsubara frequencies 
LOn. The matrix operators e, A have the structure 



/ e 
e 

-e 
VO -e 



A = 












A 








-A 








-A* 








A* 












where 



(54) 



proaches unity, i.e. monotonous increase of Ic with de- 
creasing temperature is expected in the limit of complete 
spin-flip scattering at the HS interfaces. 



2m 



(55) 



Here pfx{vfx) is the x-component of the Fermi mo- 
mentum (velocity) for a transmission channel and V{x) 
stands for the potential energy. 

Inside the half-metal it is necessary to account for 
triplet pairing which amounts to solving the two equa- 
tions 



III. GOING BEYOND QUASICLASSICS 



The quasiclassical approach used so far provides an ex- 
act solution for the problem in the limit of large number 
of conducting channels in our structure. In this case the 
quantum mechanical phases, relevant, e.g., for resonance 
effects, average out^. Such averaging is justified in a 
number of important physical situations, for instance, 
provided surface roughness of metallic interfaces exceeds 
the Fermi wavelength. On the other hand, the many 
channel limit is not the only one of physical relevance 
for the systems in question. Modern experimental tech- 
niques enable one to study the Josephson current through 
objects with few conducting channels with controllable 
change of the scattering phase^^. This renders a motiva- 
tion for calculating the Josephson current through such 
objects. In such cases the above quasiclassical formalism 
is in general insufficient, and more accurate treatment 
becomes necessary In addition, correct description of 
Andreev states in SHS junctions also requires going be- 
yond the quasiclassical approach employed in Sec. [TTl 



{iujn - e(a;i))G'ii(a:i,a;2,w„) = 6{xi - X2), (56) 

{iuJn + e{xi))G3i{xi,X2,UJn) = 0. (57) 

The current flowing through the half-metal layer is given 
by the expression 



Gii(xi,X2,UJ„ 



(58) 



Here and below an additional sum over the channel index 
k implies summation over all conducting channels of the 
junction. In the many channel limit this summation can 
be reduced to the integral over the momentum 



k 



A 



(2^) 



(59) 



Ipii \<pf 



which we already encountered in Sec. [TTl 

In what follows we will make use of the standard ap- 
proximation 



Vl[/(x)e 



] = [-PFxf{x) ± 2ipFxdxf{x)\ e 



±ipFmX 



(60) 
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which is justified for any function f{x) that varies 
smoothly on atomic distances. We further exactly fol- 
low the derivationi^. Let us fix the argument of the 
Green functions X2 inside the half-metal layer and an- 
alyze their dependence on xi. For instance, we observe 
that the solution of Eq. (f53|) decaying deep inside the 
left superconductor has the form 



Here we keep w„ > 0, denote q = exp {oJnd/2vFx) and 
define the scattering matrices 5*1.2 as non-trivial 3x3 
sub- matrices in Eq. ([?T|) . 

The solution of Eqs. ([Mll-linTl) takes the form 



(68) 



ys ^ Uizi + U2Z2, y&^Vizi + V2Z2. 



Gil 

Gil 



/2^_i ) e'P-==^ie''(^i+W2))y9(a:2). (61) 



Here yi{x2)Ty<j{x2) are arbitrary functions, a is defined 
below Eq. ([T2|) and k = ^uj'l + |Ap/w^j.. The spa- 
tial decay of the functions (G21, G31) is described analo- 
gously. 

A particular solution of Eq. ((56)) in the half-metal at 
ijJri > Q reads 



Gil ~ exp 

VFx 



WFx 



VFx 



\Xl - X2\ 



G31 = 0. 



(62) 
(63) 



Further calculation amounts to writing down the general 
solution in the half-metal layer and to matching it with 
the decaying solution in the superconducting reservoirs. 
This matching is made with the help of the scattering 
matrices, that relate outgoing and incoming waves. It is 
necessary to use two triads: The first one is composed 
of the functions Gii,G2i in the superconductor and the 
function Gn in the half- metal while the second one com- 
prises the functions G31, G41 of the superconductor and 
the function G31 of the half-metal. The second triad ac- 
counts for the hole-like excitations, hence it should be 
described by the transposed scattering matrix. 

By matching these two triads at the left and the right 
interfaces we arrive at twelve linear equations for the vari- 
ables yi,...yi2 



(64) 
(65) 
(66) 
(67) 




q ^ys 



Then the contribution to the Josephson current defined 
by the Green functions with positive Matsubara frequen- 
cies reads 



ieT J2 <l'HVi-U2) 



(69) 



w„>0,/c 



The contribution to the current from negative Matsubara 
frequencies is determined analogously. It is straight- 
forward to observe that the total Josephson current ac- 
quires the form 



/++/_= 2Rc/+, 



(70) 



i.e. it will be sufficient for our purposes to evaluate only 
the term /_|_. 

The matrices 5*1.2 relate the amplitudes of incoming 
and outgoing waves exp(±ipFxx). So if the scattering 
interface is shifted from a; = to x = a;i_2 these matrices 
are transformed as 



Sl^x 

e 



xi 

tpF^Xi Q 



Ai5i.2;=oAi, 



Ai 







(71) 
(72) 



and 



52,2: 

e 



A, 



3:2 

e-*P^ 




A252.2;=oA2 



(73) 
(74) 



B. Supercurrent 



Combining the solution of Eqs. 
we find 



(p7)l with Eqs. 
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-IX ( A{S2,a^)A*{S2,a''^)-B{S2:a^)B{S2,a^'^) \ _ ix f AjSi ,a^)A'- {Si .a''^)- BjSi ,a'^)B{Si ,a-^) ' 



/+ — ieT 



*(S2,a2)$(S2,a- 



$(Si,a2)*(Si,a-2) 



uj„>0,k e 



ix ( A{S2,a^)A'iS2,a-^)-B{S2,a^)B{S2,a-'') \ , .a'')^'- (■§! .a' 



-^)-B{Si,a^)B{Si,a-^) 
^)*(5i,a-2) 



n 



(75) 



where 



n 



a2$(5i,a2)$(52,a2) 



$(5i,a-2)$(52,a- 



(76) 



r 



and 



A{S,a^) = -S33+ 

[S31(S13S22 - S23S21) + S32(s23Sn 



S13SI2) 



+ S33(|S12P + |S21P - S11S22 - S22S*i)] + 
[■S12S2I - S11S22] [S31(S23S12 - S13S22) + 
S32(S13S21 - S23S11) + S33(silS22 - S12S21)] , (77) 



(siiS22 + S22Sn - |si2p - |s2iP) 
+ (S12S21 - S11S22) (S12S21 - 5*1^ 



(78) 



Here ip = — 2^1 — 2^2 + 2pFxd is the quantum mechanical 
phase corresponding to electron making a cycle between 
the two interfaces. As we already discussed, in the many 
channel limit it is appropriate to average the result (j80p . 
(|8ip over quickly oscillating phase ip. This averaging is 
accomplished with the aid of the relationship 



dip 



A + Be^v + 06-^"^ 2Tr - ABC' 



(82) 



It is satisfactory to observe that after such averaging in 
Eqs. ([SO)) . ([5T|) the general expression for the Josephson 
current exactly coincides with Eqs. (|151) - (|l5)) derived in 
Sec. within our quasiclassical analysis. 



^{S, a^) = S13S32 - S23S31 + 

[s3lst3('522S21 - S12S22) + S32S23(52lSn - S11S12) 

+S3iS23(|si2p - S22Sn) + S32 s J3 (sn S22 - |s2in] • 

(79) 

Making use of the parameterization ([36)) of the scattering 
matrices we eventually arrive at the general expression 
for the Josephson current 



8eT "^n sin x 



(80) 



where 



W 



{q4 [(a-i-a)(a-i 



, T , , " " , , " a cos v-\ 

X [(a^^ - a)(a"^ - acos^i^s) + |(52p] + 
q^^ [{a~^ - a){a~^ cos^ vi - a) + x 
X [(a^^ - a){a^^ cos^ 1^2 - a) + 1(52^] - 
e"^ ((a^^ - a)^ cos i^i + a\) {{a^^ - a)^ cos V2 + "2^ 



((a-i 



a 1 cos 1^1 



a) cosi/2 
-1 „^2. 



+ (5*2) ((a~^ - a)2 cos i/2 + a*^^) } 



C. Weak tunneling limit 

Let us first analyze the above general results in the 
limit of low interface transmissions i^i,;^2 ^ 1- In this 
case Eq. (|8T|) reduces to 



W 



1 



P12A4 



16 ( cosh cos p I X 

VFx 



X (c.2+A2(Reai)2)(c.2+A2(Rea2)') 



2C.2 [4 (c.2 + A2(Rea2)^) + 4 K + A^(ReaO^) 
+2.2.2(^.2 +A2)] 
+ 2iytA\lmai f {uj^ + A2(Rea2)') 

+ 4f^2 jy| A"* Im ai Re ai Im a2 Re a2 



8A2ci;2 (^iniai Re aii^i + Im a2 Re a2i/|) simysj' . 



(81) 



(83) 



The resulting expression for the current for long 
junctions is given by 
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eT 



E 

)„>0,fc 



[ul + A2(Reai)2][u;2 _^ A2(Rea2)2] (cosh. 



cos Lpk 



(84) 




FIG. 4: The function ^(99) defined in Eq. 



Here ijy{z) — dP'\TiT{z)/dz'^ is the polygamma-function. 
In the hmit T < A Eq. dHZ]) yields 



I = - 



^vlvl\l3i\\l32\vFx sin; 



We observe that at T the Josephson current grows 
with decreasing d as / cx l/d. This dependence persists 
down to d ~ Co('^i +^'2)- even smaller d <^ Co(i^i +^^1) 
we obtain 



5TA2 



uj„iyf:^^\Pi\\P2\ sinx 
(u;2+A2)3/2(z.2 + ^2)' 



(89) 



Here and below we explicitly indicate the dependence of 
the scattering phase ipk on the channel number k. At 
T = we obtain 



7C(3)e 



47rA2d3 E 

k 



i~ 121- 12 ^ '.'^fe)- 
air 



(85) 



Comparing this exact result with its quasiclassical ana- 
logue ([48]) we observe that in Eq. ((85l) the contributions 
of different channels are weighted by the function F{ipk) 
which reads (see also Fig. |4]) 



2(p 



21C(3)sin<^ 



(7r-|^|)(27r-|^|), (86) 



where — tt < (p < tt. The average of this function over 
this phase interval equals to unity. This relatively weak 
modulation of the Josephson current in our case is in 
a drastic contrast with the pronounced resonant behav- 
ior of the zero-temperature current in conventional SNS 
junctions, see, e.g., Ref. [l^. This difference is formally 
due to the presence of uj'^ in the expression for the Joseph- 
son current (|5D|) . 

Note that although in the limit d ^ ^0 resonant 
effects remain insignificant they gain importance for 
shorter junctions c? Co which will be considered be- 
low. As before, we stick to the case of weakly trans- 
mitting boundaries (small v, (3) assuming for simplicity 
Reai = Rea2 = 1- In the many-channel limit we have 
ioT i^{vl + 4) d < in 



which yields at T < A 



eA 

'2^ 



E 



vlvl\l3,\\P2\sinx 



(90) 



The above results correspond to effective averaging 
of the 93-dependent Josephson current. Let us now see 
how these results get modified if the resonance condition 
cos = \ holds for some of the conducting channels. In 
this case we have 



^^^^^E'-^ll/^i||/^2|.Lsin£ 

k 



for d > £,o{vl + V2) ^i^d 



A 

I = -eAtanh— 2^ 



j.Ml/3i|l/32|sinx 



2^2 



(91) 



(92) 



in the limit d <^ ^o{i^i + 1^2) ■ Here the sum ^'j, runs over 
the resonant channels only. 

The resonant currents (PTjl and ([^ turn out to be 
much larger than the corresponding phase-averaged con- 
tributions ([57]) . ([55)1 . Correspondingly, the resonances 
are quite narrow: Sip ^ d/^o and 6(p ^ I'l + for 
longer and shorter junctions. The off-resonant currents 
are smaller in the measure of d^/^o ^^'^ ('^i + ^2)^ 
spectively. 



eA 



327r2Td 



l + i{A/TTT) 



xy^iyfiyl\(3i\\P2\vF.sinx. (87) 



k 



D. High transmission limit 



For completeness let us analyze the case of fully trans- 
mitting interfaces, i.e. sin^ i^i = sin^ t^2 = 1- For short 
d <C Co junction the Josephson current is defined by Eq. 
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(I80ll with 



E. Andreev states 



w 



|ain«2psin^|+4^ 



"21 



A2 



(93) 



Here the initial phase is shifted by 2argai + 2arga2. 
In particular, at low temperatures T <C A and for small 
/3i^2 we have 



l/3i||/?2|sinx 



1 + I sin ^ 



(94) 



In the opposite long junction limit d ^ and at T <t; 
A the Josephson current is again given by Eq. ((80)) where 
one should substitute 



W = 



1 



IAII/32 



cosh cos ip I |q!i| \a2\ 



VFx 



2ujnd 



o, ,2 

+ ^(2+|a?| + |a2ncosh. 
A^ Wf^^ 



-2— smh (|ai| 



' A 



VFx 



(95) 



For |ai^2| ~ 1 we again recover Eq. (j85p . in which we 
should replace ai^2 by ai,2- 

Of a special interest is the case of identical interfaces 
with |/3i^2| — 1 which corresponds to equal probabilities 
for spin-up and spin-down electrons to penetrate into the 
half- metal. In this case Eqs.lUS]), (ES]) yield 



A/sgnx cos- 



for d ^ ^0 and 



7rsgnx)^Wi?^ 



(96) 



(97) 



for d » ^0- Here the Josephson phase difference is re- 
stricted within the interval — tt < x < tt and the total 
number of conducting channels in our SHS junction is 
denoted by M . One easily recognizes that Eqs. ([M)) and 
([97)1 are nothing but the 7r-shifted Kulik-Omelyanchuki^ 
and Ishii-Kulik22. current-phase relations respectively for 
short and long SNS junctions. It is also interesting to ob- 
serve that for J\f — 2 (i.e. in the single channel limit with 
two spin directions) Eq. ([97)l coincides with the result for 
the Josephson current in long SNS junctions embedded in 
a superconducting ring with odd number of electrons^. 
It also follows from Eq. ([55)1 that the current-phase rela- 
tion ((97|) for long junctions d S> sets in only in the nar- 
row parameter region |aip|Q;2p ^ Co/'^^- Correspond- 
ingly, the dependence of the Josephson current on the 
H-layer thickness changes from / cx 1 /(P to / oc 1 /d. 



The expressions for the Green functions and for the 
Josephson current derived above enable us to obtain the 
spectrum of subgap Andreev bound states inside the half- 
metal. Substituting tOn —iE we easily find the poles of 
the Green function which provide the required spectrum. 
For junctions d S> with low transparency (small /3, v) 
it is governed by the equation 



/2Ed 



cos if 



h2 cos X 



4A2(Reai)2(Rea2) 



0. (98) 



Here we assume that iJ <C A and Reai ~ Re 0:2 ~ 1. 
Though the third term in this equation is much smaller 
than the other two, it is important since it determines the 
dependence of the spectrum on the Josephson phase x- 
Note that the above equation remains valid for all values 
of if except for an immediate vicinity of the resonance 
cos ip = 1 where the full expression for W has to be taken 
into account. In the latter limit the result becomes rather 
cumbersome and is omitted here. 

It is instructive to compare the above expression with 
that for Andreev bound states in conventional SNS- 
junctions in the limit of low transmissions D12 of NS 
interfaces, see, e.g., Ref. [H. For long SNS junctions 
d ^ Co we have 



cos 



/2Ed 

V VFx 



D1D2 , , 

cos(^-l — cosx = 0. (99) 



As before, the parameter (/? is defined as 1^ = (^0 + 2pFxd. 
In the limit D12 = this equation describes particle- 
and hole- like excitations inside the normal layer with im- 
penetrable boundaries. At small but non-zero D12 these 
states get modified due to the superconducting proximity 
effect. We observe that in the presence of triplet pairing 
states inside our SHS junction Eq. ([55]) has essentially 
the same structure as Eq. ([99)) describing singlet pair- 
ing states except for an additional small factor E'^/A^ 
entering the x-dependent term in the triplet case. 

In the opposite limit of short junctions d <C Co the 
spectrum of Andreev levels can be found analogously. 
For SHS junctions we obtain 



\2 ^ \ ( m„ „ ^2 ^ 



(1 - cos^) (Reai)^ " (^6^2)^ " + 



_B2 
4A- 



:I?i2COSX = 0. (100) 



In the case Reai,2 ~ 1 this result describes excitations 
that are slightly below the gap. 

For comparison we also recall the well known expres- 
sion which defines the spectrum of Andreev levels for 
short {d <C Co) conventional SNS junctions: 



A2 



{DI + D2) 
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FIG. 5: The higher and lower Andreev levels in short SHS 
junctions with highly transparent interfaces for |/3i| = |/92| = 
0.8 



Similarly to Eq. (I100|) , this equation describes the states 
with energies slightly below the gap. We would also like 
to point out that Andreev states in SNS junctions are 
doubly degenerate. This degeneracy is in general lifted 
in the case of triplet pairing states in SHS junctions, see 
alsc^^. 

Turning to the limit of highly transparent interfaces, 
in the limit of long {d » ^o) SHS junctions for |q;i,2| 1 
we again arrive at Eq. (j98p where one should only replace 
(ReQ;i^2)^ by |ai_2p/4. In the opposite short junction 
limit d <^ £,0 we obtain the following equation, which 
is valid for arbitrary relationship between parameters a 
and (3 



4|j-(2-2|/3i||/32|cosx + |ai| 



"2 



A2" 



|aina2psin2| = 0. 



(102) 




1 V 

FIG. 6: The same as in Fig. E]for |/3i| = |/32| = 0.2 
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For positive energies < -E < A it has two solutions 
which are depicted in Figs. [5] and [6] for the case of fully 
transparent SH interfaces sin^ i/i_2 = 1. 

We observe that for \(3i,2 \ ~ 1 (see Fig. [5]) the position 
of Andreev levels significantly depends on the Josephson 
phase X whereas the ((^-dependence remains not very pro- 
nounced. In the limit |/3i,2| = 1 the energy of the lower 
Andreev level reduces to zero, while the energy for the 
upper one is determined by a simple formula 

i? = A|sin(x/2)|. (103) 

Note that this dependence is just the 7r-shifted one 
as compared to the well known dependence E — 
A|cos(x/2)| for the doubly degenerate Andreev level in 
short SNS junctions. 

In the limit of small spin-flip factors \f3i,2 \ the behavior 
of Andreev levels changes considerably, as illustrated in 
Fig. [SI In this case the dependence of the level positions 
on the Josephson phase x is much weaker, whereas their 
(^-dependence becomes more significant. 



FIG. 7: The structure of Andreev bound states in short SHS 
junctions with highly transparent interfaces for |/3i| = |/32| = 
0.8 (upper figure) and = |/32j — 0.2 (lower figure) at 
ip = 2.5 

Clearly, Andreev levels with negative energies — A < 
E < are fully symmetric with respect to the Fermi level 
and, hence, show exactly the same features, as illustrated 
by Fig. [71 We again observe that the x-dependence of all 
four Andreev states is rather pronounced for large values 
of the spin- flip parameter |/3| and it almost disappears 
for small values of \(3\. 



IV. SUMMARY 

In this paper we have developed a detailed microscopic 
theory describing dc Josephson effect and Andreev bound 
states in superconducting junctions with a half metal. 
The possibility to pass supercurrent through such SHS 
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junctions is provided by spin-flip scattering at both SH 
interfaces. As a result, superconducting correlations pen- 
etrate into the half-metal in the form of triplet pairing 
states which can carry non- vanishing supercurrent across 
the system. Our general results for the Josephson cur- 
rent, Eqs. (|43l) and ([SO]), provide a detailed description 
of the effect for sufficiently clean metals and at arbitrary 
transmissions sin^ t/i_2 of both interfaces as well as at ar- 
bitrary values of the spin-flip parameters f3i,2- In the 
tunneling limit our results reduce to those of Ref. '3 de- 
rived perturbatively in the interface transmissions. Our 
approach allows to go beyond the perturbation theory in 
1^12 and to fully account for all orders in interface trans- 
missions. This non-perturbative analysis is unavoidable 
not only at high transmissions but also in the tunnel- 
ing limit in the case of sufficiently short junctions in 
order to eliminate an intrinsic divergence of perturba- 
tive results at small d. At T = the Josephson cur- 
rent depends on the H-layer thickness as / oc for 
d » ^0 • This dependence then crosses over to / cx 1/d for 
^oii^f + i^i) < d < ^0 (see Eq. Finally, the current 

becomes ^-independent at c? ^ S.oi'^i + ^2) (see Eq. (|89| ). 
This saturation of the d-dependence of the Josephson cur- 
rent is accompanied by the change of its scaling with v 
from I IX- v'^ to I <x- . For short junctions d ^ ^0 
with few conducting channels resonant effects play an 
important role causing substantial enhancement of the 
Josephson current, see Eas. |9T|) and ([92]) . 

Our analysis demonstrates that the behavior of both 
the Josephson current and Andreev bound states in SHS 
structures crucially depends on the spin-flip parameters 
/3i,2. Similarly to Refs. 03 we observe that for incom- 
plete spin-flip scattering the temperature dependence of 



the critical Josephson current is non-monotonous with a 
maximum at non-zero T, see Fig. 3. However, this fea- 
ture disappears completely for /3i.2 1 and the Joseph- 
son current monotonously increases with decreasing T 
in this limit. Another striking feature of our results is 
that SHS junctions are characterized by the the sinu- 
soidal current-phase relation at any temperature and in- 
terface transmissions provided the spin-flip parameters 
/3i_2 are not very close to unity. Substantial deviations of 
the current-phase relation occur only provided (i) inter- 
face transmissions remain high, (ii) temperature remains 
low and (iii) the spin-flip parameters /3i,2 are very close 
to one, see Fig. 2. 

It also follows from our analysis that for symmetric 
SHS junctions (meaning that electron scattering at both 
interfaces is described by identical S'-matrices) the tt- 
junction state is usually realized. However, in a general 
case the Josephson current in SHS junctions does not nec- 
essarily show the TT-junction behavior, the current-phase 
relation is characterized by an arbitrary phase shift, see 
Eq. (|45p . It is also interesting to observe that in the 
case of fully transmitting interfaces, at T = and for 
|/3i.2| = 1 our results reduce to very simple dependencies, 
Eqs. (|96|) and ([97)) . which are essentially the 7r-shifted 
Kulik-Omelyanchuk— and Ishii-Kulik^ current-phase re- 
lations respectively for short and long SNS junctions. 

We believe that our predictions can help to identify 
triplet pairing current states in future experiments with 
SHS structures. 
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